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Abstract 

It is found that the ideal magnetohydrodjTiamic equihbrium of an axisymmet- 
ric gravitating magnetically confined plasma with incompressible flows is governed 
by a second-order elliptic differential equation for the poloidal magnetic flux func- 
tion containing five flux functions coupled with a Poisson equation for the gravi- 
tation potential, and an algebraic relation for the pressure. This set of equations 
is amenable to analytic solutions. As an application, the magnetic-dipole static 
axisymmetric equilibria with vanishing poloidal plasma currents derived recently 
by Krasheninnikov, Catto, and Hazeltine [Phys. Rev. Lett. 82 , 2689 (1999)] are 
extended to plasmas with finite poloidal currents, subject to gravitating forces 
from a massive body (a star or black hole) and inertial forces due to incompress- 
ible sheared flows. Explicit solutions are obtained in two regimes: (a) in the 
low-energy regime /5o ~ 7o ~ ~ ^ 1; where Pq, 70, Sq, and eo are related 
to the thermal, poloidal-current, flow and gravitating energies normalized to the 
poloidal-magnetic-field energy, respectively, and (b) in the high-energy regime 
/3o ~ 7o ~ (5o ~ Co ^ 1. It turns out that in the high-energy regime all four 
forces, pressure- gradient, toroidal-magnetic-field, inertial, and gravitating con- 
tribute equally to the formation of magnetic surfaces very extended and localized 
about the symmetry plane such that the resulting equilibria resemble the accre- 
tion disks in astrophysics. 



^Permanent address: University of loannina, Association Euratom - Hellenic Republic, 
Physics Department, Theory Division, GR 451 10 loannina, Greece 
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1. Introduction 



The difficult problem of equilibrium with flow has been the subject of an in- 
creasing number of investigations (e.g., in relation to the present work see Refs. 
[|I|]-|]14|, [0, pO|) both laboratory and astrophysical plasmas. Even for vanish- 



ing gravity the ideal magnetohydrodynamic (MHD) equilibrium of a symmetric 
(two dimensional) plasma with arbitrary compressible flows associated, e.g. with 
isothermal magnetic surfaces, satisfies a second-order partial differential equation 
for the poloidal magnetic flux function ip coupled with an algebraic Bernoulli 
equation for the density |^ . Depending on the value of the poloidal Mach num- 
ber (defined in Sec. 2), the above mentioned equation can be either elliptic 
or hyperbolic. Experimental evidence 0, [l^, |13|, |14|, however, conffims that (a) 
both density and temperature are to a very good approximation flux functions, i.e. 
functions constant on magnetic surfaces, and (b) the poloidal flows (e.g. involved 
in the transition from the low to the high confinement regime in tokamaks) lie 
within the ffist elliptic region. Similar conditions may also prevail in astrophysi- 
cal plasmas at least at distances not close to centers of gravity. For magnetically 
confined plasmas the equation is expected eUiptic because hyperbolicity, associ- 
ated with shock waves, would imply open magnetic surfaces and thereby abrupt 
confinement degradation. In this respect, equilibria with incompressible flows 
for which the differential equation becomes always elliptic has been of particular 
interest for both laboratory 0, |^, ^, |^, |^ and astrophysical H, ^ |10] plasmas. 

In previous studies |^ we considered ideal MHD equilibria with incom- 
pressible flows in cylindrical and toroidal geometries. For a toroidal plasma we 
found that the differential equation decouples from the pressure relation, thus 
making the problem analytically solvable. Several classes of analytic solutions of 
linearized versions of the above mentioned differential equation associated with 
sheared flows were also obtained. For vanishing flows this equation reduces to 
the Grad-Schliiter-Shafranov equation. Analytic solutions of the nonlinear Grad- 
Schliiter-Shafranov equation for a plasma with vanishing poloidal current at either 
low or high pressure confined by a dipolar magnetic field were obtained recently 
in Ref. |]15|. These studies were then extended to equilibria with purely toroidal 
flow |jl6| (being inherently incompressible because of axisymmetry), to gravitat- 
ing magnetic dipolar plasmas without flow |T^, and to plasmas with anisotropic 
pressure [jl8|- 
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The purpose of the present work is twofold: (a) to extend our equihbrium 
equations to the case of a gravitating plasma with incompressible flows, and 
(b) with employment of the separable eigenvalue technique introduced in Refs 
[|l^-[|18] to derive analytic magnetic dipolar equilibria for a plasma at finite pres- 
sure and poloidal current with incompressible sheared flows having non-vanishing 
toroidal and poloidal components, under the exertion of gravitational forces from 
a massive body (a star or a black hole). 

In Section 2 the equilibrium equations for an axisymmetric gravitating mag- 
netically confined plasma with incompressible flows are derived. These equa- 
tions are then reduced further for a plasma confined by the magnetic field of a 
point dipole in Section 3. Analytic magnetic dipolar solutions are constructed 
in Sections 4 and 5 in the following regimes: (a) in the low-energy regime 
/3o ~ 7o ~ 5o ~ eo ^ 1, where /3o, 7, ^o, and eo are related to the thermal, poloidal- 
current, flow, and gravitating energies normalized to the poloidal- magnetic-field 
energy, respectively, and (b) in the high-energy regime /5o ~ 7o ~ ~ ^ 1- 
Finally, the conclusions are summarized in Sec. 6. 

2. Stationary equilibrium equations for a gravitating 

plasma 

The ideal MHD equilibrium states of a gravitating plasma with flow are gov- 
erned by the following set of equations, written in convenient units: 



V ■ (pv) = (1) 

p(v V)v=j X B- VP-pVfi (2) 

V'^] = Gopt (3) 

V X E = (4) 
VxB=j (5) 

V ■ B = (6) 
E + V X B = 0. (7) 



Here, f2, AirGo, and pt are, respectively, the gravitation potential, the constant 
of gravity and the total density including contributions from the plasma itself 
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and from external mass sources. Standard notations are used in the rest of Eqs. 
(i-®- 

For an axially symmetric plasma the divergence-free fields can be expressed 
in terms of the functions I{R, z), F{R, z) and 0(-R, z) as 



and 



B = /V0 + V0 X V^, 
j = A>V0 - V0 X V/ 

pv = 0V0 + V0 X VF. 



(9) 



(10) 



Here, i?, 0, z are cylindrical coordinates with z corresponding to the axis of sym- 
metry, constant ip surfaces are the magnetic surfaces, and A* = i?^V ■ (V/i?^). 

Eqs. (|l|)-(0) can be reduced by means of certain integrals of the system, which 
are shown to be flux functions. To identify three of these quantities, the time 
independent electric field is expressed by E = — V$ and the Ohm's law (|^) is 
projected along V0, B and Vip, respectively, yielding F = F{ip), <l> = and 



[IF' -Q) = 



A fourth flux function is derived from the component of the force-balance equation 
(D along V0: 



/ 1- 



P 



+ WF'^' = X(V'), 



(12) 



the flux function X{ip) being related to the toroidal magnetic field (see also Eq. 
(21) in Sec. 3). Note that the toroidal quantities I{ip,R) and Q{iIj,R) are not 
flux functions. With the aid of the flux functions identifled the components of 
Eq. (0) along B and perpendicular to a magnetic surface are put in the respective 
forms 



B 



^ V Q^ VP 

V2 P / p 







(13) 



and 



V- 



p 



i?2 



+ 



F"F' 



P 



+ 



2 ^{Q/pr 



(14) 
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It is pointed out that Eqs. ( p!3[ ) and ( |I^ ) hold for any equation of state. 

We now consider incompressible flows, V ■ v = 0, which on account of Eq. (|l]) 
implies that the density is a flux function. With the aid of p = p{ip), Eq. ([TB|) 
can be integrated to yield an expression for the pressure, i.e. 

P = P.(,,)_,^_ + n__^j, (15) 

where Psi'ip) is a flux function which for vanishing flow and gravity (v = i7 = 0) 
represents the static pressure. Defining the Alfven velocity associated with the 
poloidal magnetic field, v\p = \Vip\'^/p, and the Mach number = Vp/v\p = 
[F'Y / p and inserting the expression ([T5|) into Eq. ([T^), we arrive at the elliptic 
differential equation 

(1-M2)A>--(M2)'|V^|2 

Eq. ([T6|) contains the arbitrary flux functions F{ip), ^{ip), X{ip), p{jp) and Psijp) 
which must be found from other physical considerations. For = Eqs. ( ]T5| ) and 
(|T6|) constitute concise forms of the equilibrium equations we derived formerly |p 
[Eqs. (19) and (22) therein]. We note here that the physically plausible class 
of equilibria with isothermal magnetic surfaces, T = T(?/'), was examined in our 
previous studies Q for vanishing gravity. Incompressible flows and T = T{%1)) 
imply that the pressure becomes a flux function. For toroidal plasmas it was 
shown that equilibria of this kind are possible. In addition, steady states with 
incompressible flows and isobaric magnetic surfaces were proposed in Ref. [3] as 
approximate equilibria for the Joint-European- Torus tokamak. It is also noted 
that for non-ideal plasmas with arbitrary flows, i.e. when a finite resistivity is 
introduced in Eq. (7), some of the integrals found in this section in the form of 
flux functions are destroyed and the tractability of an extension of the present 
investigation becomes questionable. 

Recapitulating, the axisymmetric ideal MHD equilibria of a gravitating plasma 
with incompressible flows are governed by the set of Eqs. (§) and (|16]) for ■?/' and 
n, and relation (|TB|) for the pressure. 
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Under the transformation |T^, |20 



(17) 



Eq. ([T6| ) reduces (after dividing by (1 — M^)^/^) to 




Eq. (|T^) is free of the nonhnear term 1/2(M^)'| V'?/'^ and, therefore, for < 1 
the equihbrium can be determined from the more tractable set of Eqs. (^, ([15|) 
and (|ll). 

3. Magnetic-dipole equilibrium equations 

The equihbrium of a plasma confined by the magnetic field of a current ring 
lying on the symmetry plane and centered at the origin of the system of co- 
ordinates is now considered. Employing spherical coordinates r, 9 and with 
/X = cos6' and R = rsin6', we seek separable solutions of Eq. ([18D of the form 



U{r,fi) = UoH{fi) 



(19) 



Here, H is an unknown function of fi alone such that H{0) = 1, and Uo and 
tq are normalization constants specifying a reference flux surface location. The 
parameter a plays the role of an eigenvalue of equations (|^) and (p!8| ) . It equals 
unity or —2 in the vacuum limit to recover the dipolar solutions ipyac oc (1 — /i^)/r 
and ipvac (1 — yU^)r^ describing, respectively, the flux surfaces far away from and 
close to the origin. We are interested in configurations symmetric with respect 
to the symmetry plane. Accordingly, the boundary conditions 



H{n — > 1) cx 1 — /i. 



dH 

djji 







(20) 



are chosen to keep the magnetic field finite at = and parallel to the axis of 
symmetry at ^ = 7r/2, respectively. Using Eqs. (g), (|n|), ([T2|), ([l^), and we 
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find for tlie magnetic field associated with Eq. (|I9] 
1 



B 



r(l-/i2)i/2 



+Bo C-^T^ (1 - M^)-V2 [(1 - ^r'^'Hi^)ee + 

V r / a ajj, 



(21) 



where Bq = aUo/vQ. Note the finite toroidal magnetic field induced by the flow 
and the poloidal currents as can also be seen from Eq. (|l^). 

We now consider a plasma subject only to gravitating forces from a star or 
black hole of mass Mg placed at r = 0; the plasma self gravity is neglected. 
Consequently, Eq. (|^) decouples from (0) and has the solution 



n 



Airr 



(22) 



Inspection of Eq. ([TsD with the gravitation potential (|2^ ) implies that the sepa- 
rable solution ([T^) is only possible provided 



Ps - PsO (^^ 



2+4 /a 



M2 



U 



jj s 2+2/a 



2+3/a 



and 



(23) 
(24) 
(25) 

(26) 



where Pso^ -^o^ Po; and $o are normalization constants associated with the refer- 
ence flux surface Uq. Inserting Eqs. (pJ]) and (|2^)- (|26|) into Eq. ([18|) , we obtain 



d^H a{a + 1) 
H : —n 



-(3oa{2 + a)H^+^/' 



;i + a)(l-/i')-^i/^+2/° 



-6oa{3 + a)(l - ^')i/^+^/" - eo«(3/2 + a)H^+^'''. 

(27) 
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Here, f3o, 70, 5o, and eo are related, respectively, to the static thermal, poloidal- 
current, flow and gravitating energies normalized to the poloidal- magnetic-field 
energy on the reference surface: 

0^3^ _ {Xo/rof/2 M%ro/Uol/2 _ poGoM,/(47rro) 

(28) 

Solutions of Eq. (^) will be constructed in the low-energy regime, /3o ~ 70 ~ 
5q ^ eo <^ 1, and in the high-energy regime, /3o ~ 70 ~ 5o ~ ^ 1- 

4. Solution in the low-energy regime 

(^0 « 70 - ^0 ~ eo <Cl) 



For this case it is convenient to put Eq. ( pTf ) in the form 



d_ 



2\2 d H 



djji \1 — ijfi 



(1 - a)(2 + a)H = -(3oa{2 + a)(l - fi^)H^+^/'' 

3 

(29) 



-7o«(l + «)i/^+'/" - 5o«(3 + «)(!- /i')'if^+^/" - eo(^ + a)(l - ^?)H^^'" 



where — 1 — //^ as /5o ~ 70 ~ (5o ~ — ^ 0. With the use of the boundary 
conditions ( pUD , integration of Eq. ( PP| ) from /i = to /i = 1 yields 

(2 + a) [(1 - a)Pi - a/3oP2] - a 7o(l + oc)P^ + 5o(3 + a)?^ + eo (^^ + a) P5 = 

(30) 

with 

Jo Jo Jo 

P4 = /x')'i/'+'/" rf/i, P5 = /'(I - /i2)i^'+'/" f^/i. 

JO JO 

(31) 

To appreciate the impact of finite pressure, finite poloidal current, fiow, and 
gravity on the vacuum equilibrium, the relations H = 1 — fi^ and a — 1 are 
employed into Eq. (^0|) except for the term 1 — a. The departure of a from unity 
is then found 

512 ^ 16 131072, 320 , , 

1-a = 3o H 7o H H eo- (32) 

1001^ 35 ^ 230945 " 693 ^ ^ 
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Therefore, the modifications of the vacuum equihbrium from the finite pressure, 
poloidal current, flow, and gravity are of the order of magnitude of /5o ~ 70 ~ 
5o ~ eo- As a resuh, an analytic solution of Eq. ( PP]) can be derived by using the 
vacuum solution H = 1— /i^ in the five terms in which H appears undifferentiated. 
Using the boundary condition ( pO|) at fi = 1 and introducing t = 1 — /i^ = sin^ 9, 
integration of Eq. (|29|) from 1 to /i yields 



di VT 



4t2(i_t)i/2 



3(1 -a) 



xdx 



(l-x)l/2 



3/3o 



* dx 

(l-x)l/2 



/■* X dx X dx 5 r 



/o (1-X)V2 



(33) 



Evaluating the integrals in Eq. (^) and integrating again, using = 0) = 1, 
we obtain the following low-energy solution valid at all distances from a point 
dipole: 



H 



1 - 



192 

Tool' 



160 
2002' 



20 








(1 


429 




572^ 


■ 6 








1(1- 


.35 


70^ 



130 



:i -t^ 



/5o 



7o 



49152 



+ 



230945 
2688 



t) 



8192 



104975 
40 



92378 
192 



™Ki-t3) + i^a-t- 



138567 



923780 



9690 



^(l-t') + — (l-t«) 
226r ^152^ ' 



5o 



100 2n 25 

(1 -t) + (1 - t^) + (1 - t^) 

23V ^ 1386^ ^ 594^ ' 



176 



Co- 



(34) 



For vanishing poloidal current, flow and gravity (70 = <5o = = 0) Eq. 
reduces to the low-pressure static equilibrium solution of Ref. [jl5| [Eq. (12 
therein] . 

5. Solution in the high energy regime 

(/3o « 70 « ^0 ~ eo > 1) 
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From the low-energy solution we anticipate that for /3o ~ 7o ~ (^o ~ eo ^ 1 it 
holds that |a| ^ 1 and, consequently, assume that 

< |a| < 1, (35) 







to be verified a posteriori. The term a{l + a)H/{l — fi^) in Eq. (pTj) is small 
everywhere, via relations (^) and (|35D , and can be neglected. 

We first search for solutions in the region 0</i</ic = l//3o'=^l- In this 
region the approximations 1 — /i"^^!,!!^! can be made and, therefore, Eq. 
can be written in the form 

^ -a f 2/5o + 70 + 3(5o + ^to) . (36) 



Integrating twice Eq. ( P^D with dH/dfi\^^Q = and i^(;U = 0) = 1 we obtain the 
solution 

H = l + a(2f3o + io + 35o + ^eo) f 1 - ^ V (37) 



Evaluation of Eq. { p7\ ) at /i = yields 

|a| = 0(l-iJ(/i,)) < 1, (38) 

consistent with assumption (PH). 

We now consider the whole regime of variation of /i. For < fi < 1 the 
terms on the RHS of Eq. (|27|) are large at /i = and rapidly decrease to zero 
as H decreases from H{fi = 0) = 1 toward H{fi = 1) = since |a| <^ 1. In 
particular, the variations of the terms 7o(l — jj^)^^a{l + a)H^^'^^" and So{l — 
/i^)Q;(3 + a)if^+^/° from 1 — /i^ are much weaker than those from H^^'^l°' and 
respectively. Therefore, in the above terms, 1 — /i^ can be approximated 
by unity and, consequently, Eq. (pT]) can be written in the form 

= -a (2/5oif^+^/" + 70^'+'/" + 36oH'+'/^ + ^eoH'+'/"^ . (39) 

Multiplying Eq. (|39D by dH/dfi and integrating from /i = where dH/djji = 0, 
we find 



-eo 



/5o f 1 - H'+'/") + 70 f 1 - + <5o f 1 - 



1 - H^+'/^)V^' + 0{J\^\. (40) 
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Integration again from fi = 0, where H{fi = 0) = 1, to /i yields 

.2+4/«^ + 70 (l 







\a\jj = 


/[ 

Jh l 



X 



X 

2+3/a 



2+2/q 



-1/2 



dx 



+6o (l - + eo (l - 

^(l-i/)(/5o + 7o + 5o + 6o)-^/^ 
To satisfy if (/i = 1) = 0, Eq. (|4l|) requires 

|«| = (/So + 70 + 5o + eo)"'/' + C(l/(/5o + 7o + ^o + eo))- 



(41) 



(42) 



Note that for /^o ~ 7o ~ 5o ~ eo S> 1, ii = 1 — /U holds everywhere except in a 
small region < ;U < (Pq + 7o + <5o + eo)^'/' ^ 1, where H remains close to unity, 
but with a large second derivative of the order of (/?o + 7o + ^o + eo)'/' [see Eq. 
(|39D]. Eq. ( ^2]) implies that the distance between adjacent flux surfaces at the 
symmetry plane /i = increases as either of Pq, ■Jq, 6o, and eo increases. Indeed, 
as a decreases the spacing must adjust to keep U ~ (ro/r)" fixed and, therefore, 
the magnetic surfaces become more extended and localized about the symmetry 
plane. The resulting equilibria resemble the accretion disks in astrophysics. It is 
noted that Krasheninnikov and Catto |l^, [l^ came to the conclusion that in the 
strong gravity limit (eo ^ 1), gravity and flow affect the flux surfaces but not 
the eigenvalue a. This does not contradict our conclusion because it concerns 
different density, current and flow regimes, i.e. arbitrary density profiles, vanish- 
ing poloidal currents, and purely toroidal flows were considered in Refs. |T^, |T^ 
while the present study concerns finite poloidal currents and magnetic surfaces of 
constant density associated with incompressible flows with non vanishing toroidal 
and poloidal components. 

6. Conclusions 

It has been shown that the equilibrium of a gravitating axisymmetric mag- 
netically confined plasma with incompressible flows is governed by a second-order 
elliptic differential equation for the poloidal magnetic flux function [Eq. (plGf)] 
containing five flux functions coupled with a Poisson equation for the gravita- 



tion potential, and an algebraic relation for the pressure [Eq. ([T5|)]. The above 
mentioned elliptic equation can be transformed to one [Eq. (|18D ] possessing a 
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differential part identical to that of the Grad-Schliiter-Shafranov equation, which 
permits the derivation of analytic solutions. 

Analytic solutions for a plasma confined by a dipolar magnetic field and sub- 
ject to gravitating forces from a massive body have been obtained in two energy 
regimes: (a) in the low-energy regime /3o ~ 7o ~ 5o ~ ^ 1; where /3o, 7o, 
6q, and eo are related to the thermal, poloidal-current, flow and gravitating en- 
ergies normalized to the poloidal-magnetic-field energy, respectively, and (b) in 
the high-energy regime /3o ~ 7o ~ ~ co ^ 1- These solutions generalize the 
static magnetic-dipole equilibria with vanishing poloidal currents obtained in Ref. 
[ p!5[| . It turns out that in the high-energy regime all four forces, pressure- gradient, 
toroidal-magnetic-field, inertial and gravitating, contribute equally to the forma- 
tion of magnetic surfaces very extended and localized about the symmetry plane 
such that the resulting equilibria resemble the accretion disks in astrophysics. 

Finally, it may be noted that, in addition to their astrophysical concern, the 
equilibrium investigations of Refs. [^, [T^, [T^ and of the present work may help 
in developing possible novel magnetic confinement devices. In this view further 
studies on the impact of compressible flows or /and self gravity on the equilibrium 
properties of plasmas confined in dipolar magnetic fields are of particular interest. 
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